Axion Cosmology 



Pierre Sikivie 



o 
o 



Theoretical Physics Division, CERN, CH-1211 Geneve 23, Switzerland 
and 

■ Department of Physics, University of Florida, Gainesville, FL 32611, USA 
| sikivie@phys.ufl.edu 

o 

^ | ■ Outline: 

1. Thermal axions 

2. Axion field evolution 

3. The domain wall problem 
CS| ■ 4. Cold axions 

5. Axion miniclusters 

■ 6. Axion isocurvature perturbations 

o 



Note: For background information on the Strong CP Problem, its resolution 
| by an axion, and the laboratory and astrophysical constraints on axions, the 

reader is referred to the lecture notes by Roberto Peccei and Georg Raffclt in 
this series. Here, we will be concerned only with the cosmological properties 
of axions. Different authors may use different definitions of the axion decay 
constant. In the present notes the axion decay constant is represented by the 
^ ■ symbol f a and defined by the action density for QCD plus an axion 

£ Q CD+a = -\G%G^ V + l -d^a 

+ J2 e(*7 M ^ - m q )q +^(6 + ^)G%G a ^ (1) 

>5Z7r J a 

where a is the axion field before mixing with the r\ and it mesons. Eq. ([1]) uses 
standard notation for the chromomagnetic field strengths, the strong coupling 
constant and the quark fields. The axion mass (after mixing with the rj and 
7r mesons) is given in terms of f a by 

f lO 12 GeV \ 

m a ~ 6 [ieV I J . (2) 

f a is related to the magnitude v a of the vacuum expectation value that breaks 
the t/pQ(l) symmetry: f a = v a /N. N is an integer characterizing the color 
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anomaly of J7pq(1). N = 6 in the original Peccei-Quinn- Weinberg- Wilczek 
axion model. All axion couplings are inversely proportional to /„. 



1 Thermal axions 

Axions are created and annihilated during interactions among particles in the 
primordial soup. Let us call "thermal axions" the population of axions estab- 
lished as a result of such processes, to distinguish them from the population 
of "cold axions" which we discuss later. 

The number density ?i* (i) of thermal axions solves the Boltzmann equa- 
tion (Jl|): 

iin th 

-g- + 3ff< h = r(rC - < h ) (3) 

where 

r = ^2 n t <<j lV > (4) 

i 

is the rate at which axions are created and annihilated, Hit) is the Hubble 
expansion rate, and 

rC = (5) 

is the number density of axions at thermal equilibrium. £(3) = 1.202.. is the 
Riemann zeta function of argument 3. In Eq. ((J), the sum is over processes 
of the type a + i <-> 1 + 2, where 1 and 2 are other particle species, rii is 
the number density of particle species i, cr^ is the corresponding cross-section, 
and the brackets indicate averaging over the momentum distributions of the 
particles involved. 

Unless unusual events are taking place, T cx R~ 1 where R(t) is the scale 
factor, and Eq. ([5]) implies therefore 

dn cq 

+ 3Hn? = . (6) 
Combining Eqs. ijHJ) and one obtains 

|[i? 3 « h - rC)] = -r R a (nf < q ) ■ (7) 

Eq. |(7J) implies that a thermal distribution of axions is approached exponen- 
tially fast whenever the condition 

r > h (8) 

is satisfied. So, we have a thermal population of axions today provided in- 
equality © prevailed for a few expansion times at some point in the early 
universe, and provided the thermal population of axions thus established did 
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not subsequently get diluted away by inflation, or some other cause of huge 
entropy release. 

The least model-dependent processes for thermalizing axions are: 1) a + 
q{q) 9 + 2) a + g <-> q + q and 3) a + g «-> g + g. The corresponding 
diagrams are shown in Fig. 1. These processes involve only the coupling of the 
axion to gluons, present in any axion model [see Eq. {TJ], and the coupling 
of quarks to gluons. A detailed treatment is given in ref. ([2). We give only a 



a) b) 



P 




c) 



Fig. 1. Processes which produce thermal axions in the early universe 
rough estimate here. The processes of Fig. 1 have cross-sections of order 

2 

where a s = At temperatures T > 1 TeV, the densities of quarks, anti- 
quarks and gluons are given by 

n q = n- q = 27 ^ T 3 

„ = 16 ^ T 3 . (10) 
The Hubble rate is given by the Friedmann equation 

# 2 - — P = — + gA/}(T)) — T , (11) 

where A/i (T) (A//(T)) is the total effective number of bosonic (fermionic) spin 
degrees of freedom at temperature T. For T > 1 TeV, 

N = N b + -Af f = 107.75 (12) 



4 Pierre Sikivic 



if we assume no new degrees of freedom other than those of the Standard 
Model plus an axion. Combining everything and setting a c ~ 0.03, one finds 

\^ r T i GeV\ 2 T 

g * H ~ 2 HH H)^GeV ■ (13) 

Thus we find that the processes of Fig. 1 keep axions in thermal equilibrium 
with the primordial soup till the temperature 

T °> ~ 5 ' 1011GCV (W^ev) 2 ■ < 14 > 

Note that the calculation is not valid when T > v a = N f a , since the PQ 
symmetry is restored then. In particular, in view of Eq. (|14|) . the processes 
under consideration only produce a population of thermal axions if /„ < 
2 N 10 12 GeV. 

We will see in Section 4 that f a has to be less than approximately 10 12 GeV 
to avoid overdosing the universe with cold axions. That bound suggests that 
the processes of Fig. 1 do produce a population of thermal axions. We should 
keep in mind however that this thermal axion population may be wiped out by 
a period of inflation with reheat temperature less than Tq s . So it is interesting 
to search for processes that may re-establish a thermal axion population later 
on. We briefly discuss two such processes. 

First we consider the Compton-like scattering process Q(Q) + g <-> Q(Q) + 
a. The quark Q may be a known quark or a new heavy quark. There are re- 
lated processes in which the gluon is replaced by a photon or Z-boson, and/or 
the quark is replaced by a lepton. However Q(Q) + g «-> Q(Q) +a is in some 
sense the least model dependent among the Compton-like processes since, in 
every axion model, there is at least one colored fermion Q which carries PQ 
charge, and hence to which the axion couples. The cross-section is 

o-Q ~ a s (^^j 7^2 for T > m Q 

~ a s \ for T < ma . (15) 

The relevant regime is when T > raq since the Q number density is Boltzmann 
suppressed for T < tuq. Using a c — 0.05 , N = 107.75, we have 

rr i m Q /!0 12 GeV\ 2 

n Q < a Q v > H~ l ~ r(2 ,J GeV) {—^- ) (16) 

for T > mQ. So this process produces a population of thermal axions provided: 

2 



1 ^ 2 ' 107GeV (ll)^) ^ 
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The axions decouple then at a temperature Tpc ~ m Q- 

Let us also consider the process tt + tt <-> 7r + a. Since the axion 
necessarily mixes with the 7r°, this process is model-independent as well. It 
occurs at T ~ 200 MeV, after the QCD phase-transition but before the pions 
have annihilated. The cross-section is of order 



1 

72 

./ a 



(18) 



Using TV = 17.25, we find 



TT , /3- 10 8 GeV\ 
< ovw > ir 1 - ( - J (19) 

at T ~ . The 7r + 7r <-> 7r + a process has the advantage of occuring very 
late, so that any thermal axion population it establishes cannot be wiped out 
by inflation. (Inflation occuring that late would also wipe out the baryons.) 
However, Eq. (fT9|) indicates that it is ineffective unless the bound f a > 10 9 
GeV from SN1987a is saturated. 

We have seen that, under a broad set of circumstances, a population of 
relic thermal axions is produced. For f a > 10 9 GeV, the axion lifetime exceeds 
by many orders of magnitude the age of the universe. Between their last 
decoupling, at temperature TYj, and today the thermal axion population is 
merely diluted and redshifted by the expansion of the universe. Their present 
number density is 

,_ C(3)^ 3 (Rn 



"«e°> = ^ r Hi£J (20) 

where ^ is the ratio of scale factors between the time tjj of last decoupling 
and today. Their average momentum is: 

(pf(t )) = — ^-——Te)—^- = 2.701 To—— . (21) 

If (p^ (to)) >> ni a , the energy distribution is thermal with temperature 

T a0 = T D ^ . (22) 

If there is no inflation, nor any other form of entropy release, from to till the 
present, T a o is related to the present cosmic microwave background temper- 
ature T 7 o = 2.735 K by the conservation of entropy. Taking account of the 
fact that electron-positron annihilation occurs after neutrino decoupling, one 
finds 

The average momentum of relic thermal axions is 
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i 

(^(t ))=2.110- 4 cv(^)" , (24) 



and their number density is 

+„\ - 

cm 3 V Af D 



th/ . 7.5 / 106.75 \ 
»- (*«>) - — s • (25) 



2 Axion field evolution 

The thermal axions discussed in the previous Section are quantum fluctua- 
tions about the average background value of the axion field. The evolution of 
the average axion field, from the moment Upq(1) gets spontaneously broken 
during the PQ phase transition to the moment the axion acquires mass during 
the QCD phase transition, is the topic of this Section. 

The Upq(1) symmetry gets spontaneously broken at a critical temperature 
Tpq ~ v a , where v a is the vacuum expectation value of a complex field <j>{x). 
The action density for this field is of the form 

C^ l -d^d^-- A (^4>-vlf + .... (26) 

where the dots represent interactions with other fields in the theory. At T > 
Tpq, the free energy has its minimum at <fi — 0. At T < Tpq, the minimum is 
a circle, whose radius quickly approaches v a as T decreases. Afterwards 

(<f>(x)) = v a e'-W/"" (27) 

where a(x) is the axion field before mixing with the 7r° and rj mesons. a(x) 
has random initial conditions. In particular, at two points outside each other's 
causal horizon the values of a are completely uncorrelated. 

It is well-known that the size of the causal horizon is hugely modified dur- 
ing cosmological inflation. Without inflation, the size of the causal horizon 
is of order the age t of the universe. But, during an inflationary epoch, the 
causal horizon grows exponentially fast and becomes enormous compared to 
t. There are two cases to consider. Case 1: inflation occurs with reheat tem- 
perature smaller than Tpq, and the axion field is homogenized over enormous 
distances. The subsequent evolution of this zero momentum mode is relatively 
simple. Case 2: inflation occurs with reheat temperature larger than Tpq. In 
case 2, in addition to the zero mode, the axion field has non zero modes, and 
carries strings and domain walls as topological defects. 

The early universe is assumed to be homogeneous and isotropic. Its cur- 
vature is negligible. The space-time metric can therefore be written in the 
Robertson- Walker form: 



- ds 2 = dt 2 - R(t) 2 dx ■ dx 



(28) 
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where x are co-moving spatial coordinates and R(t) is the scale factor. The 
equation of motion for a(x) in this space-time is: 

D^a(x) + VMx)) = (d 2 + 3^d t - ^l)a(x) + Vfa(x)) = , (29) 

where V a is the effective potential for the axion field, and prime indicates 
a derivative with respect to a. V a results from non-perturbative QCD effects 
associated with instantons (0). They break Upq(1) symmetry down to a Z(N) 
discrete subgroup (0). V a {a) is therefore periodic with period Aa — 2irv a /N = 
2irf a . We may write such a potential qualitatively as 

V a = f a m 2 a (t)[l-co S (^)] , (30) 

J a 

where the axion mass m a (t) — m a (T(t)) is a function of temperature and 
hence of time. Eq. ((27)) implies that the axion field has range a € [0, 2irv a }- 
Hence there are N degenerate vacua. The discrete degeneracy implies the 
existence of domain walls, which will be discussed in Section 3. 

Substituting Eq. (|30| into Eq. (|29|) . the equation of motion becomes 

(d? + 3|a, - -^Vl)a(x) + ml(t)f a sin(^) = . (31) 

The non-perturbative QCD effects associated with instantons have amplitudes 
proportional to 

e -^m ~ (^)U-I^ (32 ) 



where Nf is the number of quark flavors with mass less than T. Eq. (f32|) 
implies that the axion mass is strongly suppressed at temperatures large com- 
pare to the QCD scale, but turns on rather abruptly when the temperature 
approaches v1q CD 

Because the first three terms in Eq. (f3~Tj) are proportional to t~ 2 , the axion 
mass is unimportant in the evolution of the axion field until m a (t) becomes 
of order j. Let us define a critical time t\. 

m„(ti)ti = 1 . (33) 

The axion mass effectively turns on at t\. m a (T) was obtained (0;@; 
a calculation of the effects of QCD instantons at high temperature 
result is: 

,„„ (T) . 4 . 10 - 8e v(i^X)(^) 4 ,34, 

when T is near 1 GeV. The relation between T and t follows from Eq. (fTTj) 
and H — i . The total effective number M of thermal spin degrees of freedom 
is changing near 1 GeV temperature from a value near 60, valid above the 
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quark-hadron phase transition, to a value of order 30 below that transition. 
Using Af ~ 60, one has 

,20 1 ( t V /10 12 GeV N 



which implies: 



-.(0 = 0.7 10"-^) {—f— ) , (35) 



The corresponding temperature is: 



Tl ~ i GeV I ^ ) . (37) 

Eq. (|3"S")) implies -4ln(m a (£)) < m a (t) after ti. So, at least for a short while 
below 1 GeV, as long as Eq. ([34]) remains valid, the axion mass changes 
adiabatically. The number of axions is the adiabatic invariant. Conservation 
of the number of axions after t\ allows us to estimate the energy density 
of axions today from an estimate of their number density at t\. When the 
temperature drops well below 1 GeV, the dilute instanton gas calculations 
which yield Eq. (fM)) are no longer reliable. Complicated things happen, such 
as the confinement and chiral symmetry breaking phase transitions. However, 
because m a >> H then, it is reasonable to expect the number of axions to be 
conserved, at least in order of magnitude. 



2.1 Zero mode 

In case 1, where inflation occurs after the PQ phase transition, the axion field 
is homogenized over enormous distances. Eq. (f3Tj) becomes (0; S 0) 

^ + Is)«w + '»"<^ 8i "(x)^ ■ < 38 » 

where we used R(t) cx \ft. For t << t\, we may neglect m a - The solution is 
then 

a(t) = a + ait~^ (39) 



where clq and ai are constants. Eq. (|39[) implies that the expansion of the 
universe slows the axion field down to a constant value. 

When t approaches ii, the axion field starts oscillating in response to 
the turn on of the axion mass. We will assume that the initial value of a is 
sufficiently small that f a sm(a/f a ) ~ a. Let's define tp: 

a(t) = t-§ip(t) . (40) 



The equation for ip(t) is 
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^+^(t))m = o (4i) 

where 

^(t)=ml(t) + ^ . (42) 

For t > ti, we have 4? lna> « U) m a . That regime is characterized by the 
adiabatic invariant ipQ(t)u>(t), where ipo(t) is the changing oscillation ampli- 
tude of ip(t). We have therefore 

m ~ ^£= cos ( f dt'Lo(t')) , (43) 

where C is a constant. Hence 

a{t)=A{t)cos(j dt'oj(t')J , (44) 

with 

^) = ™ • ( 45 ) 
y/m a (t) U 

Hence, during the adiabatic regime, 

A 2 (t)m a (t) rxt'i (x R(t)~ 3 . (46) 

The zero momentum mode of the axion field has energy density p a — ^m^A 2 , 
and describes a coherent state of axions at rest with number density n a = 
■^m a A 2 . Eq. (|46|) states therefore that the number of zero momentum axions 
per co-moving volume is conserved. The result holds as long as the changes 
in the axion mass are adiabatic. 

We estimate the number density of axions in the zero momentum mode 
at late times by saying that the axion field has a random initial value a(ti) = 
f a oti and evolves according to Eqs. (|44l l45f for t > ti. cci is called the 'initial 
misalignment angle'. Since the effective potential for a is periodic with period 
27i" fa j the relevant range of a\ values is —it to +7r. The number density of zero 
momentum axions at time t\ is then (0; 0; 0) 

nT c - (h) ~ im a (t 1 )(a(t 1 )) 2 = ^K) 2 (47) 

where we used Eq. (|3"3"|) . We will use Eq. (|47p in Section 4 to estimate the zero 
mode contribution to the cosmological energy density of cold axions. 

A more precise treatment would solve Eq. (|38p for t ~ t\, e.g. by nu- 
merical integration, to obtain the exact interpolation between the sudden 
(t < fi) and adiabatic (f > ti) regimes. An additional improvement is to solve 
Eq. ([38| without linearizing the sine function, thus allowing large values of a\. 
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Although these improvements are desirable, they would still leave the num- 
ber of axions unknown in case 1 because the initial misalignment angle a.\ is 
unknown. In case 2, the zero mode contribution to the axion number density 
is also given by Eq. (|47p but the misalignment angle ot\ varies randomly from 
one horizon to the next. 



2.2 Non zero modes 

In case 2, where there is no inflation after the PQ phase transition, the axion 
field is spatially varying. Axion strings are present as topological defects, and 
non zero momentum modes of the axion field are excited. In this subsection, 
we consider a region of the universe which happens to be free of strings. Strings 
will be added in the next subsection. 

The axion field satisfies Eq. (f3Tj) . We neglect the axion mass till t ~ t\. 
The solution of Eq. (f3Tj) is a linear superposition of eigenmodes with definite 
co-moving wavevector k: 

<z(x,*) = J d 3 k a(k,i)e jk x (48) 

where the a(k,i) satisfy: 

5 t 2 + ^ t + ^a(k,i) = 0. (49) 

Eqs. ([28]) and (j48]) imply that the wavelength X(t) = ^-R(t) of each mode 
is stretched by the Hubble expansion. There are two qualitatively different 
regimes in the evolution of a mode, depending on whether its wavelength is 
outside (A(i) > t) or inside (A(i) < t) the horizon. 

For X(t) ^> t, only the first two terms in Eq. (|49|) are important and the 
most general solution is: 

a(k, t) = a (k) + a_ 1/2 (k)t- 1 / 2 . (50) 

Thus, for wavelengths larger than the horizon, each mode goes to a constant; 
the axion field is so-called "frozen by causality" . 

For X(t) < t, let a(k,t) = Rr§ (i)^(k, t). Eq. becomes 

(a t 2 + ^ 2 (t))^(k,i) = o (5i) 

where 

^ )= i^) + i-i?%- (52) 

Since ^ ln(w) ui, this regime is again characterized by the adiabatic invari- 
ant ipoO 1 , t)w(t), where ipo(k, t) is the oscillation amplitude of ip(k,t). Hence 
the most general solution is: 
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a(k,t) = -^-rcosf / dt'u(t')) , (53) 



R(t) 

where C is a constant. The energy density and the number density behave 
respectively as p a ,u ~ %zjjj oc and n 0ik - ^p a ,k oc indicating that 

the number of axions in each mode is conserved. This is as expected because 
the expansion of the universe is adiabatic for modes with X(t)t <C 1. 

Let us call ^^(uj,t) the number density, in physical and frequency space, 
of axions with wavelength A = 2£ for uj > t _1 . The axion number density in 
physical space is thus: 

n„(i)= / dtu ^(w,t) , (54) 
J t -i dw 

whereas the axion energy density is: 

f dfi 

Pa(t) = / dui—^(u>,t)u. (55) 
J t -i dui 

Under the Hubble expansion axion energies redshift according to ui' = ui-jjr, 

and volume elements expand according to AV' = AV (j^j j whereas the 
number of axions is conserved mode by mode. Hence 

dn a ( R'\ 2 dn a R 

d7 ( "' <) = UJ ^r (w ^'°- (56) 

Moreover, the size of for uj ~ j is determined in order of magnitude by 
the fact that the axion field typically varies by v a = Nf a from one horizon to 
the next. Thus: 



dria t ±\ A 
uj— — (uj, t)Auo 

dui 



dn a (\ \ fly I, N9 lN 2 tf , s 



From Eqs. (J5BJ) and ([S7]). and i? oc -ft, we have (0) 

dn Q/ s TV 2 / 2 , x 

d^-ad- (58) 

Eq. (|58|) holds until the moment the axion acquires mass during the QCD 
phase transition. 



2.3 Strings 

In case 2 axion strings are present as topological defects in the axion field 
from the PQ to the QCD phase transitions (|l(J). The energy per unit length 
of an axion string is 
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H = nv 2 a ln(v L) . (59) 

L is an infra-red cutoff, which in practice equals the distance to the nearest 
neighbor string. Because they are strongly coupled to the axion field, the 
strings decay very efficiently into axions. We will see that practically all axions 
produced by string decay are non-relativistic after t\ . Because each such axion 
contributes m a to the present energy density, it is important to evaluate the 
number density of axions emitted in string decay. This is our main goal in this 
subsection. 

At a given time t, there is at least on the order of one string per horizon. 
Indeed the axion field is completely uncorrelated over distances larger than t. 
Hence there is non-zero probability that the random values of a(x, t) wander 
from zero to 2irv a along a closed path in physical space if that closed path 
has size larger than t. When this is the case, a string perforates the surface 
subtended by the closed path. 

At first, the strings are stuck in the primordial plasma and are stretched 
by the Hubble expansion. During that time, because R(t) oc y/i, the density 
of strings grows to be much larger than one per horizon. However expan- 
sion dilutes the plasma and at some point the strings become unstuck. The 



temperature at which strings start to move freely is of order (jlll ) 



T * ~ 2 107GeV (lOlfeev) ■ (60) 

Below T*, there is a network of axion strings moving at relativistic speeds. 
Axions are radiated very efficiently by collapsing string loops and by oscil- 
lating wiggles on long strings. By definition, long strings stretch across the 
horizon. They move and intersect one another. When strings intersect, there is 



a high probability of reconnection, i.e. of rerouting of the topological flux (|12f). 
Because of such 'intercommuting', long strings produce loops which then col- 
lapse freely. In view of this efficient decay mechanism, the average density of 
long strings is expected to be of order the minimum consistent with causality, 
namely one long string per horizon. Hence the energy density in long strings: 

PstT (t)=^^^Apl Hvat) , (6i) 

where £ is a parameter of order one. 

The equations governing the number density n s ^ T (t) of axions radiated by 
axion strings are (jlll ) 

^£ = -2Hp stI - ( 6 2) 
dt dt v ' 



and 



-dT-- 311 ^ (63) 
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where u>(t) is defined by: 

1 1 f dk dp s 



At) rf P°fa->° J k dt dk 

w dt J 



(64) 



k is axion momentum magnitude. dps d " t " a (t) is the rate at which energy den- 
sity gets converted from strings to axions at time t, and ffi *7l a (t, k) is the 
spectrum of the axions produced. ui(t) is therefore the average energy of 
axions radiated in string decay processes at time t. The term —2Hp s tr — 
+Hp str — 3-f/pstr in Eq. ([62]) takes account of the fact that the Hubble expan- 
sion both stretches (+Hp str ) and dilutes (—3Hp str ) long strings. Integrating 
Eqs. (|6"1"1- |6"5)) . setting H = Jr, and neglecting terms of order one versus terms 
of order ln(v a t), one obtains 

r , (66) 



t2 Jt PQ t'*Lj{t>) 

where ipQ is the time of the PQ transition. 

To obtain n s ^(t) we need to know ui(t), the average energy of axions 
radiated at time t. If u)(t) is large, the number of radiated axions is small, and 
vice-versa. Axions are radiated by wiggles on long strings and by collapsing 
string loops. Consider a process which starts at ti n and ends at t nn , and which 
converts an amount of energy E from string to axions. t- ln and ta n are both 



taken to be of order t. It is useful to define the quantity (jl3f ) 



N ax (t) = J dk^(t)l . (66) 

where k is wavevector, and ^§(t) is the wavevector spectrum of the a field. 
At the start (t = tin), only string constributes to the integral in Eq. (|6"6")) . 
At the end (t = tan), only axions contribute. In between, both axions and 
string contribute. The number of axions radiated is N a = -/V ax (£fi n ), and their 
average energy is lu = -M-. The energy stored in string has spectrum ^ oc \ 
for k min < k < fc max where fc max is of order v a and fc m ; n of order 2zc ^ St. If 
I = is the length of string converted to axions, we have 

AUfrn) = utv \ , • (67) 



Hence 



w ln(u t)fc min 

where r is the relative change in -/V ax (i) during the process in question: 



(68) 



— /V (+• ") ' 1 j 

J v ax V^in J 
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fc m ; n is of order ^£ where L is the loop size in the case of collapsing loops, 
and the wiggle wavelength in the case of bent strings. L is at most of order t 
but may be substantially smaller than that if the string network has a lot of 
small scale structure. To parametrize our ignorance in this matter, we define 
X such that the suitably averaged fc m ; n = X^f- Combining Eqs. (|65|) and (|68|) 
we find: 

PrN 2 f 2 

n -tr( t )~$ Lsl , (70 ) 

x t 

where f is the weighted average of r over the various processes that convert 
string to axions. One can show |28l ) that the population of axions that were 
radiated between ipQ and t have spectrum cx -p- for j < k < —j= — , 

irrespective of the shape of -^fff , provided t » tpQ. 

At time t\, each string becomes the edge of N domain walls, and the 
process of axion radiation by strings stops. Since their momenta are of order 
t^ 1 at time t\, the axions radiated by strings become non-relativistic soon 
after they acquire mass. We discuss the string decay contribution to the axion 
energy density in Section 4. For the moment we turn our attention to the 
domain walls which appear at t\ in case 2. 



3 The domain wall problem 

Axion models have an exact, spontaneously broken, discrete Z(N) symme- 
try. Z(N) is the subgroup of Upq(1) which does not get broken by non- 
perturbative QCD effects Q). The spontaneously broken Z(N) symmetry 
implies a N fold degeneracy of the vacuum. The N vacua are at equidis- 
tant points on the circle at the bottom of the 'Mexican hat' potential for the 
Peccei-Quinn field </). An axion domain wall is the minimum energy field con- 
figuration which interpolates between neighboring vacua. Note that there are 
axion domain walls even when N = 1. In this case both sides of the domain 
wall are in the same vacuum (indeed there is only one vacuum) but the in- 
terpolating field configuration winds around the bottom of the Mexican hat 
potential once. The properties of walls in N = 1 models are for most purposes 
identical to those of walls in N > 2 models. The only seemingly important 
difference is that the walls of N = 1 models are quantum mechanically unsta- 
ble. Even so, their decay rate per unit surface and time is exponentially small 
and negligible in the context of our discussion. 

When the axion mass turns on, at time t\, each axion string becomes the 
edge of N domain walls. The domain walls produce a cosmological disaster 
unless there is inflation after the PQ phase transition (case 1), or unless N = 1. 
Indeed, let's consider the implications of case 2 if N > 2. Since there are two or 
more exactly degenerate vacua and they have identical properties, the vacua 
chosen at points outside each other's causal horizon are independent of one 
another. Hence there is at least on the order of one domain wall per causal 
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horizon at any given time. In case 2, the size of the causal horizon is of order 
t, the age of the universe. Thus the energy density in domain walls 

PwW > j (71) 

where a is the wall energy per unit surface, given by (ji; 14) 

a ~ 9f 2 a m a ~ 5.5 • 10 10 GeV 3 ( -J^—\ . (72) 



IQ 12 GeV, 

The energy density in axion domain walls today (to — 14 Gyr) 



MWaj-o.io-^^^j (73) 

would exceed by many orders of magnitude the critical energy density, of order 
10 -29 gr/cm 3 , for closing the universe. This would be grossly inconsistent with 
observation. Let's see what would happen. 

Let t w be the age of the universe when the domain walls start to dominate 
the energy density. The condition H 2 ~ ^pPw and Eq. (jTT) imply 

Domain walls are gravitationally repulsive (jlEl ; [1(1 [13) ■ They accelerate away 
from each other with acceleration 2irGa and, after a time of order (27rG<7) , 
recede at the speed of light. By averaging over volumes containing many cells 
separated by walls, the equation of state of a wall dominated universe is seen 
to be 

2 

Pw = -gPw • (75) 

Conservation of energy 

d( Pw R 3 ) = ~ Pw d{R 3 ) , (76) 

where R is the scale factor, then implies p w oc i. This scaling law and the 
Friedmann equation 

H = [r) = — p - (77) 

imply that a domain wall dominated univers expands according to 

R cx t 2 . (78) 

The domain wall dominated universe has an accelerated expansion. One may 
be tempted to attribute the present acceleration of the expansion of the uni- 
verse (3; 3) to domain walls. However, a domain wall dominated universe 
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is far less homogenous than ours. It would be divided into cells separated by 
rapidly expanding walls. Inside each cell, concentrated near the cell's center, 
would be a clump of matter and radiation with total energy of order 

M ~ P {t w )tl ~ 1O U M (^) . (79) 

One cannot identify these clumps with galaxies because neighboring clumps 
are receding from one another at close the speed of light. 

There are three solutions to the axion domain wall problem. The first 
solution is to have inflation with reheating temperature less than the PQ 
phase transition temperature, i.e. postulate case 1. In case 1 the axion field is 
homogenized by inflation and there are no strings or domain walls, and hence 
no domain wall problem. The second solution is to postulate N = 1. The third 
solution is to postulate a small explicit breaking of the Z(N) symmetry. The 
viability of the second and third solutions is less obvious. We discuss them in 
succession. 

3.1 N = 1 

The above arguments, showing the existence of a domain wall problem, arc 
valid only when the vacuum is multiply degenerate. They do not apply to the 
N = 1 case. On the other hand, since N = 1 models contain domain walls too, 
it is not immediately clear that they are free of difficulties. However, N = 1 
is a solution (fiol : [ill ) as we now discuss. 

In the circumstances under consideration (case 2) axion strings are present 
in the early universe from the time of the PQ phase transition to that of the 
QCD phase transition. At temperature T\ each string becomes the boundary 
of a single domain wall. To see what the network of walls bounded by string 
looks like, a cross-section of a finite but statistically significant volume of the 
universe near time t\ was simulated in refs. (j2ll : l9l). The simulation shows 
that there are no infinite domain walls which are not cut up by any string. 
The reason for this is easily understood. An extended domain wall has some 
probability to be cut up by string in each successive horizon it traverses. The 
probability that no string is encountered after traveling a distance I along the 
wall decreases exponentially with I. 

The question now is: what happens to the network of walls bounded by 
strings. The walls are transparent to the thermalized particles in the primor- 
dial soup, whose typical momentum is of order one GeV, but the walls have 
a large reflection coefficient for non-relativistic axions such as the cold axions 
that were produced by vacuum realignment and wall decay (0). The drag 
on the motion of the walls which results from their reflecting cold axions is 
important at time t\ but turns off soon afterwards as the cold axions are di- 
luted by the expansion of the universe ([9]). The wall energy per unit surface 
a(t) ~ 8m a (t)fa is time dependent. [Eq. ([72]) is for zero temperature.] The 
string at the boundary of a wall is embedded into the wall. Hence its infra-red 
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cutoff L, in the sense of Eq. ([55]) . is of order the wall thickness m a 1 Q). The 
energy per unit length of such string is therefore 

/.~7r/ a 2 ln(/ a /m a ) . (80) 

The surface energy E a of a typical (size ~ t\) piece of wall bounded by string 
is <j(t)t\ whereas the energy in the boundary is ~ /iij. There is a critical 
time t% when the ratio 

E a (t) ^ 8m a (t)ti 
it^ 7rln(/ a /m a ) 

is of order one. Using Eqs. (|35p and (|36p. one estimates: 

t2 - 10 " sec (^4v) 1/3 (82) 

Zl0 12 GeV\ 1/6 
T 2 ~ 600MeV f — - — j . (83) 

After £2 the dynamics of the walls bounded by string is dominated by the 
energy in the walls whereas before t 2 it is dominated by the energy in the 
strings. A string attached to a wall is pulled by the wall's tension. For a 
straight string and flat wall, the acceleration is: 

er(i) 8m a (t) m a {t) 1 m a (t) 

<M*J - - ; 77 7 7 — oq ~ 7 7TT- 

fi 7rln(/ a /m a ) 23 tim a (t 2 ) 

Therefore, after t 2 , each string typically accelerates to relativistic speeds, in 
the direction of the wall to which it is attached, in less than a Hubble time. 
The string will then unzip the wall, releasing the stored energy in the form of 
barely relativistic axions. We estimate in section 4.3 how much walls bounded 
by string contribute to the present cosmological axion energy density. 

A very small portion of the domain wall energy density is in walls which 
are not bounded by string, and which form closed surfaces such as spheres 
or torii |2~0: lit). Such closed walls do not decay by the process just described. 
Instead, they oscillate and emit gravitational waves (fiol). Using the quadrupole 
formula, we may estimate the gravitational wave power emitted by a closed 
wall of size I oscillating with frequency to ~ £~ 1 : 

p „ „ G(o .^)V „ Ga 2 f (85) 

at 

Eq. ([55)1 implies the lifetime: 

T gra v ~ (Ga)- 1 ~ 2 ■ 10 3 sec ^ 1 ° 12 / GCV ^) , (86) 



independently of size. 
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3.2 A small breaking of PQ symmetry 

A third solution to the domain wall problem is to postulate a small explicit 
breaking of the Z{N) symmetry, and hence of PQ symmetry (0). The symme- 
try breaking must lift completely the degeneracy of the vacuum and be large 
enough that the unique true vacuum takes over before the walls dominate 
the energy density. On the other hand, it must be small enough that the PQ 
mechanism still works. This solution does not appear very attactive and we 
will see below that there is little room in parameter space for it to occur, but 
it is a logical possibility. 

To get rid of the walls, we add to the RHS of Eq. (26| a tiny U?q(1) 
breaking term which lifts the vacuum degeneracy completely, e.g.: 

5V = -£((f>e- zS + h.c.) . (87) 

To add such a term to an axion model by hand seems rather unnatural. 
However it is conceivable that a small Z7pq(1) breaking term is in fact a 
natural property of the ultimate theory. This would be the case, for example, 
if the low energy effective theory at some energy scale has an automatic PQ 
symmetry which is broken in the full theory. Be that as it may, Eq. ([57]) yields 
a small correction to the effective potential for the axion field: 

6V a = -2v a i cos (j- - 6 J . (88) 

The unique true vacuum is the one for which | 8 — -p | is smallest. Its energy 
density is lowered by an amount of order £v relative to the other now quasi- 
vacua. As a result, the walls at the boundary of a region in the true vacuum 
are subjected to an outward pressure of order £,v a . Since the walls are typically 
a distance t apart, the volume energy £ v a i 3 associated with the lifting of the 
vacuum degeneracy grows more rapidly than the energy at 2 in the walls. At 
a time r ~ |2- , the pressure favoring the true vacuum starts to dominate the 
wall dynamics and the true vacuum takes over, i.e. the walls disappear. The 
energy stored in the network of walls and strings decays into gravitational 
waves (0). The true vacuum must take over before the walls dominate the 
energy density. Using Eq. (|74p . we obtain: 

r^'^^f^C). (89) 

& a \ fa J 

On the other hand £ is bounded from above by the requirement that 5V does 
not upset the PQ mechanism. SV shifts the minimum of the effective potential 
for the axion field, inducing a 9 ~ . The requirement that 9 < 10~ 10 
implies: 

10 sec f fa \ fnn , 
T ~N \ lO^GeV J ■ (90) 



Eqs. (|89|) and ([90]) indicate that there is little room in parameter space for 
this third solution to the axion domain wall problem, but it is not ruled out. 
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4 Cold axions 

Our main purpose in this Section is to estimate the cosmological energy den- 
sity in cold axions, and estimate their velocity dispersion. The discussion in 
the previous two sections identified the following sources of cold axions: 

- vacuum realignment 

- zero momentum mode 

- higher momentum modes 

- axion string decay 

- axion domain wall decay. 

In case 1, only the contribution from vacuum realignment with zero momen- 
tum is present. In case 2, all sources are present. 

4.1 Vacuum realignment 

Zero momentum mode 

Let us start with case 1 which is easiest to analyze. The axion number 
density at time t\ is given by Eq. (|47|) in terms of the initial vacuum misalign- 
ment angle ot\. In case 1, where the universe was homogenized by inflation 
after the PQ phase transition, ol\ has the same value everywhere. We saw in 
Section 2 that, barring any sudden changes in the axion mass during the chi- 
ral symmetry breaking phase transition, the number of axions is an adiabatic 
invariant after t\. Hence the axion density from the vacuum realignment zero 
mode at a later time t 

< ac '°W~^-K) 2 (§-) 3 for easel (91) 

where ^ is the ratio of scale factors between t\ and t. Thus we find 

P T'°{h) ~ ^(«i) 2 (I;)' ^ case 1. (92) 

for the axion energy density today. 

In case 2, the initial misalignment angle ct\ is different from one QCD 
horizon to the next. Since the average of (ai) 2 over many QCD horizons is of 
order one, we have 

pr^o)~^(f£) 3 for case 2. (93) 

However in case 2, there are additional contributions. 



Higher momentum modes 
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In case 2, the axion field is not constant even within each horizon volume. 
It has wiggles inherited from earlier epochs when the horizon was smaller 
and the axion field was inhomogeneous on correspondingly shorter scales. 
The associated density of axions in physical and frequency space is given in 
Eq. (|58p . Integrating over uj > , we find the contribution from vacuum 
realignment involving higher momentum modes 

N 2 f 2 

< ac,1 (ti)~— ±± for case 2. (94) 

2ti 

Almost all these axions are non-relativistic after t\. Hence 

pjr ,i (to) „ rn^n forcase2 . (95) 

Note that, except for the factor TV 2 , p™ ' 1 and p™ c '° are °f the same order of 
magnitude. 



4.2 String decay 

The evolution of axion strings between the PQ and QCD phase transitions was 
discussed in section 2.3 . The number density of axions emitted in string decay 
is given by Eq. (|70[) in terms of quantities £, f , and x, which we introduced 
to parametrize our ignorance of various aspects of string evolution and decay. 
Because their spectrum oc p- over the range j < k < —j= — , the bulk 



/MPQ 

of axions emitted in string decay have momenta of order t^ 1 at time t\, and 
become non-relativistic soon after they acquire mass. Therefore, the string 
decay contribution to the axion energy density today is 

pf (to) = manfih) (|i) 3 - rnj-^^l ^ . (96) 



We now discuss the factors on the RHS of Eq. (|96|) which are specific to the 
string decay contribution. 

£ This parameter determines the density of the string network, Eq. (|61[) . 
with £ = 1 corresponding to a density of one long string per horizon. In ref. 
(fl3h it was argued that £ ~ 1 because global strings can decay efficiently into 
axions, and therefore the number density of long strings should be close to the 
minimum consistent with causality. In numerical simulations of global string 



networks in an expanding universe (|22f). it was found that indeed £ ~ 1. So 



there appears to be good ground for setting £ ~ 1. 

X This parameter defines the low wavevector edge of the ^g- oc -p- 
spectrum through fc m i„ = X~r- X an d £ are related since the average interstring 
distance controls both. On dimensional grounds, x ^ \/£- So the effect of small 
scale structure in the axion string network partially cancels out in the RHS 
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of Eq. x is expected to be of order one, but the uncertainty on this is at 
least a factor two. 

f This parameter defines the average energy of the axions emitted in 
string decay, through Eqs. (|68l69p . It is the unknown on which most of the 
debate has focused in the past. Two basic scenarios have been put forth, which 
we call A and B. The question is: what is the spectrum of axions radiated by 
strings? The main source is closed loops of size L ~ t. Scenario A postulates 
that a bent string or closed loop oscillates many times, with period of order 
L, before it has released its excess energy and that the spectrum of radiated 
axions is concentrated near =£. In that case one has f ~ ln(v a ti) ~ 67. 
Scenario B postulates that the bent string or closed loop releases its excess 
energy very quickly and that the spectrum of radiated axions is 4f °^ \ with a 
high frequency cutoff of order 2irv a and a low frequency cutoff of order . In 
scenario B, the initial and final spectra ^§ of the energy stored in the axion 
field are qualitatively the same and hence f ~ 1. In scenario A, the string 
decay contribution dominates over the vacuum realignment contribution by 
the factor ln(v ii), whereas in scenario B the contributions from string decay 
and vacuum realignment have the same order of magnitude. 

Many authors ((H E3) have argued in favor of scenario A, 
adopting the point of view that global strings are similar to local strings 
and that their coup ling to the axion field can be treated perturbatively. My 
collaborators and I (| 1 ll; 1 1 3f) have argued in support of scenario B, emphasizing 
that the dynamics of global strings is dominated by the energy stored in the 
axion field and that there is no reason to believe that this energy would 
behave in the same way as the energy stored in the string core. The numerical 
simulations of the motion and decay of axion strings in refs. (fl3t [28h give 
strong support to scenario B. These simulations are of oscillating strings with 
ends held fixed, of collapsing circular loops, and of collapsing non-circular 
closed loops with angular momentum. Over the range of ln(v a L) accessible 
with present technology [2.5 < hi(v a L) < 5.0], it was found that r ~ 0.8 
for closed loops and r ~ 1.07 for oscillating strings with ends held fixed. No 
dependence of r on \n(v a L) was found for closed loops, and for bent strings 
with ends held fixed r was found to slightly decrease with increasing ln(v a L), 
whereas scenario A predicts r to be proportional to h\(v a L). 



4.3 Wall decay 

The final contribution to the cold axion cosmological energy density in case 2 is 
from the decay into non-relativistic axions of axion walls bounded by string. 
We assume here that N = 1. Indeed if TV > 1, the domain wall problem 
is presumably solved by introducing a small breaking of PQ symmetry, as 
described in section 3.2. In that case, the axion walls decay predominantly 
into gravitational radiation (@). 

Let ts be the time when the decay effectively takes place and 7 = m ^ t3 ) 
the average Lorentz 7 factor of the axions produced, to' is their average energy. 
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The density of walls at time ti was estimated to be of order 0.7 per horizon 
volume (Q). Hence the average energy density in walls is 

«.„.(«) - 0.7 ^ (f ) 3 - (0.7)(8)m a (i)| j^f (97) 

between t\ and t%. We assumed that the energy in walls simply scales as a(t). 
After time £3, the number density of axions produced in the decay of walls 
bounded by strings is of order 

d.w. / 1\ Pd.w. 



Note that the dependence on £3 drops out of our estimate of n^' w '. In the 
simulations of the motion and decay of walls bounded by string in ref. ([91) it 
was found that 7 ~ 7 for ln(-^-) ~ 4.6 but that 7 increases approximately 
linearly with hi(y/\v a /m a ). If this behaviour is extrapolated all the way to 
]n(y/Xv a /m a ) — 60, which is the value in axion models of interest, then 7 ~ 60. 
In that case the contribution from wall decay is subdominant relative to those 
from vacuum realignment and string decay. 

4.4 The cold axion cosmological energy density 

To estimate the cosmological energy density of cold axions in case 2, we neglect 
the contribution from wall decay and assume that scenario B is correct for 
the string contribution. By adding the RHS of Eqs. ([9"3]l. ([95]) . and (|96j) with 
N = r = £ = x — 1) we nn d 

Pa(*o) ~ 2 ^ f m a for case 2. (99) 

Eq. (j9"2"|) gives the cold axion cosmological energy density in case 1. To de- 
termine the ratio of scale factors R\/ we assume conservation of entropy 
from time t\ till the present. The number M\ of effective thermal degrees of 
freedom at time t\ is of order 60. Keeping in mind that neutrinos decouple 
before electron- positron annihilation, one finds 

!) 3 ,„. 06 3(^)\ (100, 
Combining Eqs. ([36]). (f99]) . ([92]) and flOO]) . and dividing by the critical density 

^~°- 15 (i^cv) (t) al forcasel 
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where h is defined as usual by H = h 100km/s-Mpc. 

Eqs. (|10ip are subject to many sources of uncertainty, aside from the un- 
certainty about the contribution from string decay. The axion energy density 
may be diluted by the entropy release from heavy particles which decouple 
before the QCD epoch but decay afterwards |29t [301 [3lh. or by the entropy 
release associated with a first order QCD phase transition. On the other hand, 
if the QCD phase transition is first order (jH; [H; [3; El) , an abrupt change of 
the axion mass at the transition may increase J? a ■ A model has been put forth 
|36|) in which the axion decay constant f a is time-dependent, the value f a (ti) 
during the QCD phase-transition being much smaller than the value f a today. 
This yields a suppression of the axion cosmological energy density by a factor 
( /yO )2 conl p are d to the usual case. Finally, it has been proposed that the 
axion density is diluted by 'coherent deexcitation', i.e. adiabatic level crossing 
of m a (t) with the mass of some other pseudo-Nambu-Goldstone boson which 
mixes with the axion (1371). 



4.5 Velocity dispersions 

The axions produced by vacuum realignment, string decay and wall decay 
all have today extremely small velocity dispersion. In case 1, where the ax- 
ions are produced in a zero momentum state, the velocity dispersion is zero. 
(This ignores the small quantum mechanical fluctuations created during the 
inflationary epoch, which will be discussed in Section 6.) 

In case 2, we distinguish two subpopulations of cold axions: pop. I and 
pop. II, with the second kind having velocity dispersion typically a factor 10 3 
to 10 larger than the first. The pop. I axions are those produced by vacuum 
realignment or string decay and which escaped being hit by moving domain 
walls. They have typical momentum pi(ti) ~ j- at time t\ because they are 
associated with axion field configurations which are inhomogeneous on the 
horizon scale at that time. Their velocity dispersion is of order: 

1 ( R A o in -i7 A0- 5 eV\ 5/6 i? o 



m a ti \li ) \ m a J K 

The corresponding effective temperature is of order 0.5 10~ 34 K ( 10 eV )5 
today. This very cold, indeed! 

Pop. II are axions produced in the decay of domain walls and axions that 
were hit by moving domain walls. Axions produced in the decay of domain 
walls have typical momentum ^11(^3) ~ 7^10(^3) at time £3 when the walls 
effectively decay. Their velocity dispersion is therefore of order: 

m a (t 3 )R 3 _ in _ 13 _ flO- 5 eV\ 1/6 R 



m a R \ m a J R 

where q = 7 "'° n ^ 3 " > ^7 parametrizes our ignorance of the wall decay process. 
We expect q to be of order one but with very large uncertainties. There is 
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however a lower bound on q which follows from the fact that the time t 3 when 
the walls effectively decay must be after t 2 when the energy density in walls 
starts to exceed the energy density in strings. Using Eq. (|82l) . we have 

,, 7^3) ^3 , im a {t 2 ) R 2 _ 7 / 10- 5 cV \ 2/3 
m a Ri m a i?i 130 \ m a J 

Since computer simulations suggest 7 is of order 60, pop. II axions have 
much larger velocity dispersion than pop. I axions, by a factor 10 3 or more. 
Whereas pop. II axions are relativistic or near relativistic at the end of the 
QCD phase transition, pop. I axions are definitely non-relativistic at that time 
since m a » l/t\. The axions which were produced by vacuum realignment 
or string decay but were hit by relativistically moving walls at some time 
between t\ and £3 should be included in pop. II since they are relativistic 
just after getting hit. The next Section will highlight the differences in the 
behaviours of the two populations of cold axions. 

The very low velocity dispersion of cold axions, and their extremely weak 
couplings, imply that these particles behave as cold collisionless dark matter 
(CDM). CDMparticles lie at all times on a 3-dim. hypersurface in 6-dim. 



phase-space (38; .39). As a result, CDM forms discrete flows and caustics. The 
number of discrete flows at a given physical location is the number of times the 
3-dim. hypersurface covers physical space at that location. At the boundaries 
between regions with differing number of flows, the 3-dim. hypersurface is 
tangent to velocity space. The dark matter density is very large on these 
surfaces, which are called caustics. The density diverges at the caustics in the 
limit of zero velocity dispersion. 



5 Axion miniclusters 

If there is no inflation after the PQ phase transition (case 2), the initial mis- 
alignment angle a± changes by 0(1) from one QCD time horizon to the next. 
Hence, the fluid of cold axions produced by vacuum realignment is inhomo- 
geneous with = 0(1) at the time of the QCD phase transition. As will 
be shown shortly, the streaming length of pop. I axions is too short for these 
inhomogcncities to get erased by free streaming before the time t e q of equal- 
ity between matter and radiation, when density perturbations start to grow 
in earnest by gravitational instability. At time i oq , the = 0(1) inho- 
mogeneities in the axion fluid promptly form gravitationafly bound objects, 
called axion miniclusters |40l: [4lk I42I lit) . The properties of axion miniclusters 
are of concern to experimentalists attempting the direct detection of dark 
matter axions on Earth. Indeed those experiments would become even more 
challenging (than they are already) if most of the cold axions condense into 
miniclusters and the miniclusters withstand tidal disruption afterwards. Of 
course, these issues only arise in case 2. There are no axion miniclusters in 
case 1. 
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As described above, there are two populations of cold axions, pop. I and 
pop. II, with velocity dispersions given by Eqs. (|102|) and (|103[) respectively. 
Both populations are inhomogeneous at the time of the QCD phase transition. 
The free streaming length from time t\ to t cq is: 

t { = R(t eq )J^dt^~(3(ti)^t^\n(^} . (105) 

The time of equality and the corresponding temperature are respectively t cq ~ 
2.3 10 12 sec and T eq ~ 0.77 eV. The free streaming length should be compared 
with the size 

* mc ~ii^^ v ^^2-10 13 cm(^-^) (106) 
Hi \ m a J 

of axion inhomogeneities at t eq . Using Eq. (|102|) we find for pop. I: 

^^ln(M.2.WH^y /3 . (107) 
«mc hm a \h J \ m a J 

Hence, in the axion mass range of interest, pop. I axions do not homogenize. 
At t e q most pop. I axions condense into miniclusters. The typical size of axion 
miniclusters is £ mc and their typical mass is (U 3) 



/10~ 5 eV\ 5/3 

M mc ~ r] p a (*eq)C ~ V 5 ' 1O~ 13 M , (108) 

\ m a J 

where r\ is the fraction of cold axions that are pop. I. We assumed that all pop. 
I axions condense into miniclusters, and used Eq. (| 1 1 1 — case 2) to estimate 

Pa (tcq)- 

Using Eq. (|103|) . we find for pop. II: 

£°~gln(V)~42g. (109) 



t 3 

Using Eq. (|104[) and assuming the range 7 ~ 7 to 60, suggested by the nu- 
merical simulations of ref. (0), we conclude that pop. II axions do homogenize 
and hence that the axion energy density has a smooth component at t e q- 

However, pop. II axions may get gravitationally bound to miniclusters 
later on. It seems rather difficult to model this process reliably. A discussion 
is given in ref. (|9|). It is concluded there that the accretion of pop. II axions 
results in miniclusters which have an inner core of pop. I axions with density 
of order 10 -18 gr/cm 3 and a fluffy envelope of pop. II axions with density of 
order 10 -25 gr/cm 3 . 

When a minicluster falls onto a galaxy, tidal forces are apt to destroy it. If 
a minicluster falls through the inner parts of the Galaxy (r < 10 kpc) where 
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the density is of order 10 -24 gr/cm 3 , its fluffy envelope of pop. II axions will 
likely be pulled off immediately. This is helpful for direct searches of dark 
matter axions on Earth since it implies that a smooth component of dark 
matter axions with density of order the halo density permeates us whether or 
not there is inflation after the Peccei-Quinn phase transition. Even the central 
cores of pop. I axions may eventually get destroyed. When a minicluster passes 
by an object of mass M with impact parameter b and velocity v, the internal 
energy per unit mass AE given to the minicluster by the tidal gravitational 
forces from that object is of order (lol ) 



r2«2 



C 2 M 2 

AE ~ (no) 

whereas the binding energy per unit mass of the minicluster E ~ G Pmc^mc- 
If the minicluster travels a length L — f3t through a region where objects of 
mass M have density n, the relative increase in internal energy is: 

AE Gpljt 2 . . 

— (HI) 

& Pmc 

where pu — Mn. Eq. (jllip follows from the fact that AE is dominated by the 
closest encounter and the latter has impact parameter 6 m j n : "KoL\ in nL = 0(1). 
Note that =§- is independent of M. A minicluster inner core which has spent 
most of its life in the central part of our galaxy only barely survived since 
£§■ ~ 1(T 2 in that case. 

The direct encounter of a minicluster with Earth would be quite rare, 
happening only every 10 4 years or so. The encounter would last for about 3 
days during which the local axion density would increase by a factor of order 
10 6 . 



6 Axion isocurvature perturbations 



In this Section we describe the isocurvature perturbations IH 0; H; 0; 0; 
471 ) produced if inflation occurs after the Peccei-Quinn phase transition, and 
derive the constraints on axion parameters from the absence of isocurvature 
fluctuations in CMBR observations. 

If the reheat temperature after inflation is less than the temperature Tpq 
at which C/pq(1) is restored (case 1), the axion field is present during inflation 
and is subject to quantum mechanical fluctuations, just like the inflaton. In 
fact, since the axion field is massless and weakly coupled like the inflaton, it 
has the same fluctuation spectrum (4^ 4{| 5(| 51) 
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where Hj is the expansion rate during inflation. As before, x are comoving 
spatial coordinates. The axion fluctuations described by Eq. (|112p are com- 
monly written in shorthand notation as 8a = y 1 . The fluctuation in each 
axion field mode is "frozen in" after R(t)/k exceeds the horizon length HJ 1 . 
We do not consider here the possibility of fluctuations in the axion decay 



constant f a during inflation. Such fluctuations are discussed in Refs. (|52l ; l53l ; 
0). 

At the start of the QCD phase transition, the local value of the axion field 
a(x, t) determines the local number density of cold axions produced by the 
vacuum realignment mechanism [see Eq. (|47p] : 

n„(x,ti) = A a ( Xj<1 )2 (113) 

where a(x,fi) = a(x, t\)/f a is the misalignment angle. The fluctuations in 
the axion field produce perturbations in the cold axion density 

= M£ = _Hi_ ( 

n a ax irfa,ai 

where a x = a(ti) = f a a% is the value of the axion field at the start of the QCD 
phase transition, common to our entire visible universe. These perturbations 
obey 

SpTih) = (H5) 

since the vacuum realignment mechanism converts energy stored in the quark- 
gluon plasma into axion rest mass energy. In contrast, the density perturba- 
tions produced by the fluctuations in the inflaton field 

m m ieh m Hi) 

satisfy 

= 7 ■ (H6) 

/^matter ^ Prad 

Density perturbations that satisfy Eq. are called "adiabatic", whereas 

density perturbations that do not satisfy Eq. are called "isocurvature" . 

Isocurvature perturbations, such as the density perturbations of Eq. (|115[) . 
make a different imprint on the cosmic microwave background than do adia- 
batic ones. The CMBR observations are consistent with pure adiabatic per- 
turbations. This places a constraint on axion models if the Peccei-Quinn phase 
transition occurs before inflation. 

Before we derive this constraint, two comments are in order. The first 
is that, if the Peccei-Quinn transition occurs after inflation, axion models 
still predict isocurvature perturbations but not on length scales relevant to 
CMBR observations. Indeed we saw in the previous section that in this case 
(case 2) the axion field fluctuates by order f a from one QCD horizon to the 
next. Those fluctuations produce isocurvature perturbations on the scale of 
the QCD horizon, which is much smaller than the length scales observed in 
the CMBR. Their main phcnomenological implication are the axion miniclus- 
ters which were discussed in Section 5. The second comment is that, if the 
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Peccei-Quinn phase transition occurs before inflation (case 1), the density 
perturbations in the cold axion fluid have both adiabatic and isocurvature 

components. The adiabatic perturbations (-^f- = ^ ra ^ = are produced 
by the quantum mechanical fluctuations of the inflaton field during inflation, 
whereas the isocurvature perturbations are produced by the quantum me- 
chanical fluctuations of the axion field during that same epoch. The adiabatic 
and axion isocurvature components are uncorrelated. 

The upper bound from CMBR observations and large scale structure data 
on the fraction of CDM perturbations which are isocurvature is of order 30% 
in amplitude (10% in the power spectrum) (5^; 6(1 61; ~6^, 6^; ID). Allowing 



for the possibility that only part of the cold dark matter is axions, the bound 
on isocurvature perturbations implies 

_?a_ = J^ ^_ = / a < 0.3 — , (117) 

PCDM Pa PCDM ^ ] a (X\ J/CDM Pm 

where we used Eq. (| 1 14|) . is the amplitude of the primordial spectrum 
of matter perturbations. It is related to the amplitude low multipole CMBR 



anisotropics through the Sachs- Wolfe effect (|65l ; l66t l67h . The observations 
imply ^& ~ 4.6 10" 5 (jH). 



In terms of cti, the cold axion energy density is given by Eq. (jlOip . We 
rewrite that equation here, assuming h ~ 0.7: 

n ^°- 15 (whevY al ■ (118) 



It has been remarked by many authors, starting with S.-Y. Pi (69), that it is 
possible for f a to be much larger than 10 12 GeV because a\ may be accidentally 
small in our visible universe. The requirement that S7 a < ^cdm = 0-22 implies 

^><0A( 10l \ GeV Y . (119) 



7T \ fa 

Since — tt < u\ < +tt is the a-priori range of ai values and no particular 
value is preferred over any other, |^| may be taken to be the "probability" 
that the initial misalignment angle has magnitude less than |ai|. (Strictly 
speaking, the word probability is not appropriate here since there is only one 
universe in which a\ may be measured.) If | — | = 2-10 -3 , for example, f a may 
be as large as 10 16 GeV, which is often thought to be the "grand unification 
scale" . 

The presence of isocurvature perturbations constrains the small ot\ sce- 
nario in two ways (j47h . First, it makes it impossible to have ot\ arbitrarily 
small. Using Eq. (fm]) . one can show that the fluctuations in the axion field 
cause the latter to perform a random walk (0) characterized by the property 

i / d 3 x < (Sa(x, t) - 5a(0, t)) 2 >= 4nH] ln(i?fc max ) (120) 

V Jy 
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where the integral is over a sphere of volume V = ^-R 3 centered at x = 0, 
and fc max is a cutoff on the wavevector spectrum. a\ cannot be smaller than 
the RHS of Eq. (|120p with R equal to the size of the present universe and 
fcmax equal to the Hubble rate at QCD time, redshifted down to the present. 
Since Q a < 0.22, this implies a bound on Hj. Translated to a bound on the 
scale of inflation Ai, defined by H] = ^-GAj, it is 

A,< 5- i "Gev( I5I ^)™ . (121) 

Second, one must require axion isocurvature perturbations to be consis- 
tent with CMBR observations. Combining Eqs. (| 1 1T[) and (|118|) . and setting 
J? CDM = 0.22, &a = 4.6 10~ 5 , one obtains 

yl 7 <10 13 GeV fl"* (j^^j " • (122) 

Let us keep in mind that the bounds (] 1 2 1 [) and (|122|) pertain only if the reheat 
temperature Trh < 2pq. One may, for example, have Q a — 0.22, f a ~ 10 12 
GeV, and ylj ~ 10 16 GeV, provided T RH > 10 12 GeV, which is possible if 
reheating is sufficiently efficient. 
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